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Advanced gravitational wave interferometric detectors will operate at their design sensitivity 
with nearly ~1MW of laser power stored in the arm cavities. Such large power may lead to the 
uncontrolled growth of acoustic modes in the test masses due to the transfer of optical energy to 
the mechanical modes of the arm cavity mirrors. These parametric instabilities have the potential 
of significantly compromising the detector performance and control. Here we present the design 
of “acoustic mode dampers” that use the piezoelectric effect to reduce the coupling of optical to 
mechanical energy. Experimental measurements carried on an Advanced LIGO-like test mass shown 
a 10-fold reduction in the amplitude of several mechanical modes, thus suggesting that this technique 
can greatly mitigate the impact of parametric instabilities in advanced detectors. 


I. INTRODUCTION 

The network of advanced gravitational wave detectors 
currently under construction (two LIGO [1] detectors in 
the USA, the Advanced Virgo [2] detector in Italy, and 
Kagra [3] in Japan) promises to open the new window of 
gravitational wave astronomy within this decade. 

Sensitivity to gravitational wave strains of order 
requires high optical power circulating in the arm cavities 
of these detectors. Eor instance, up to 750 kW of optical 
power will be sustained in the steady-state regime inside 
the Advanced LIGO arm cavities. 

It has been experimentally observed that the stored 
energy in a resonant cavity can leak from the optical 
modes to the mechanical modes of the cavity mirrors 
via a 3-mode interaction [4]. Given sufficiently high cir¬ 
culating optical power, and mirror materials with very 
low mechanical loss as required to avoid thermal noise, 
the uncontrolled growth of test mass acoustic modes can 
lead to Parametric Instabilities (PI) O [6] . If left unad¬ 
dressed PI will prevent high power operation, and thus 
limit the astrophysical output of gravitational wave de¬ 
tectors. While Advanced LIGO will serve as the primary 
example in this paper, all advanced gravitational wave 
detectors are susceptible to these instabilities. 

Several schemes have been proposed to damp PI [7118) . 
In particular, solutions directly applicable to Advanced 
LIGO involve active damping of acoustic modes by means 
of the test mass electro-static drive actuators [9], and 
thermal tuning of the optical modes using the test mass 
ring heaters [10]. A significant constraint on any tech¬ 
nique is that it must preserve the inherently low me¬ 
chanical loss of the test mass in the gravitational wave 
frequency band to maintain a low level of thermal noise. 

Here we present a novel method to passively control 
PI by reducing the Q-factor of the test mass acoustic 
modes with small resonant dampers. These “acoustic 
mode dampers” (AMD) dissipate the strain energy of the 
test mass mode through a resistive element after convert- 
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ing it to electrical energy via the piezoelectric effect (see 
figure . 

The resonant nature of the AMD allows it to effectively 
damp test mass acoustic modes without introducing sig¬ 
nificant mechanical loss at lower frequencies where ther¬ 
mal noise can limit detector performance. With respect 
to other proposed solutions, this approach has the ad¬ 
vantage of being simple, self-contained, and completely 
passive. Models indicate that AMDs can provide a broad¬ 
band reduction in the Q of mechanical modes relevant to 
PI, and are therefore particularly beneficial in the pres¬ 
ence of a large number of unstable modes. 

The structure of this paper is as follows. In section [TI| 
we set the stage by giving a brief overview of parametric 
instabilities, including equations of particular relevance 
to evaluating AMD performance. Section m presents a 
simple 1-dimensional model of the AMD which highlights 
the principles of AMD operation. This is followed by a 
description of the detailed finite element model (EEM) 
used to analyze the AMDs ability to suppress PI when 
attached to an Advanced LIGO test mass. The EEM 
predictions are compared with experimental results from 
a full-scale prototype in section [^ Einally, in section |V| 
we discuss an AMD design that will provide Advanced 
LIGO with protection from instability, without signifi¬ 
cantly increasing test mass thermal noise. 


II. PARAMETRIC INSTABILITIES (PI) 

The acousto-optic interactions responsible for para¬ 
metric instabilities have been extensively studied Bin 
m- They consist of a scattering process and radia¬ 
tion pressure operating together in an optical cavity in a 
closed-loop manner. The graphical representation of this 
process is shown in Eigj^ 

The e-folding growth time, or “ring-up time”, of an 
acoustic mode in the presence of a 3-mode interaction is 
r = 2Qml (<^m (7^ “ 1 )) [ 12 ] : whcrc Qm and ujm are the Q- 
factor and angular frequency of the mode, respectively, 
and R is the parametric gain. When R> 1 the amplitude 
increases exponentially until a saturation point is reached 
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FIG. 1. PI described as a positive feedback process. A 
steady state cavity field inside the interferometer arm cavity 
is distorted by scattering off a vibrating mirror surface. Two 
transverse optical sidebands are created. Both sidebands ex¬ 
ert force on the test mass via radiation pressure. When the 
energy dissipation of the acoustic mode and the rate of work 
done by the radiation pressure is unbalanced, one of the side¬ 
bands excites the exponential growth of the acoustic mode 
amplitude. The dynamic of this process is commonly de¬ 
scribed in terms of the parametric gain R, with R > 1 in 
the case of instability. 


m- The parametric gain R can be defined as 


R = An^PeAQm X X X r(Aw) (1) 

OJ^hom 01^00 


where Pc is the optical power circulating in the arm cav¬ 
ity, Uq is the optical frequency of the light, and Suqq and 
5uhom are the cavity linewidths (full width, half maxi¬ 
mum) for the fundamental and the higher-order optical 
mode, respectively. The parameter A measures the spa¬ 
tial overlap between the acoustic and the higher-order 
optical modes; F is representative of the interferometer 
optical configuration and is a function of the 3-mode in¬ 
teraction tuning Auj = ujm — 27TAuhom^ where Auhom 
is the frequency difference between the fundamental and 
higher-order optical modes. For Aoo -> 0, the parame¬ 
ter F reaches its maximum (see [14] for a more detailed 
description). 

Unstable acoustic modes with parametric gain up to 
R 100 may arise in Advanced LIGO in the 10-90 kHz 
band [5]. To prevent these instabilities, a damping mech¬ 
anism must be introduced to reduce the Q-factor of all 
unstable acoustic modes in this frequency band without 
introducing excess thermal fluctuation in the detection 
band of 10 Hz to 1 kHz. 


III. MODEL OF THE ACOUSTIC MODE 
DAMPER (AMD) 

In order to reduce the Q test mass acoustic modes we 
designed a resonant AMD which can be attached to the 
test mass and provide dissipation via the piezo-electric 
effect. 

In this section we first describe the interaction between 
the AMD and the test mass with a simple 1 -D model, 
then we present a complete Finite Element Model of the 
entire system. 


A. Simplified 1-D Model 

The AMD and test-mass system can be described as a 
pair of coupled oscillators with a large mass ratio. The 
AMD mass m is attached to the much more massive test 
mass via piezo electric shear plates, which are modeled as 
a lossy spring with complex spring constant of magnitude 
k and loss angle g. 

The test mass acoustic mode for which we would like 
to estimate the impact of the AMD is simplified in this 
model to a mass M, equal to the modal mass of the 
acoustic mode, attached to a fixed reference by a lossless 
spring K. The coupled systems is then excited by the 
radiation pressure force F applied to the TM mode, as 
shown in Fig. 

At frequencies near the resonance of the AMD, the 
lossy spring produced by the piezoelectric material and 
resistive load will dissipate the energy of the excited 
acoustic mode, as seen in Fig. 

For this system of coupled oscillators, the amplitude 
A{uj) of the acoustic mode driven by force F at angular 
frequency uj is 


A{uj) = 


^2 ^ ^2 


Muj‘^ y + fiY + + fiY 


where J = 1 — ^ e = 1 — 

o K o k ^ m 
= ^D = and/i = —. 


M 


M 


The resulting effective Q-factor is 


max(A(ci;)) 77 ^ + (1 — p)^ 
A(uj = 0) ppp 


( 2 ) 


( 3 ) 

( 4 ) 


where p = cjo/cji, and we assume /i <C 1. 

When the acoustic mode resonance is near that of the 
AMD, p ^ |1 — p|, the acoustic mode Q is reduced to 
Qeff pIp. When the acoustic mode resonance is well 
above the AMD resonance, Qeff — p/pp^ and when it is 
well below the AMD resonance, Qeff — ^IpPPi assuming 
p‘^ < 1 . 

To suppress Pis, the test mass acoustic mode Q-factors 
only need to be reduced from ^ 10^ to 10^ — 10^. Using 
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FIG. 2. Overview of the Acoustic Mode Damper (AMD). The AMD can be described as a small spring-mass system attached 
to a larger mass M. M represents a modal mass of an acoustic test mass mode which vibrates due to an external force F 
(radiation pressure). Such induced vibration causes deformation of the lossy spring of the damper. The non-zero loss angle of 
the damping spring (piezoelectric material) results in the lag of strain with respect to stress and thus energy dissipation. As a 
consequence of this dissipation process the amplitude of the acoustic mode vibration is reduced. 


this simple model we can estimate that a 3g AMD with 
f] = 0.1 on a 10kg test mass, can give Qeff ^ 10^ for 
resonances with | < p < 3. 

However, this model ignores a number of important 
factors. One of these is the location of the AMD relative 
to the nodes and anti-nodes of each test mass acoustic 
mode. Quantitatively speaking, the modal mass M of 
a given mode should be increased in this model by the 
ratio of the displacement at the AMD location to that 
of the mode’s anti-node squared M' = M(xmax/^AMD)^- 
Thus, an AMD located near a node will have a reduced 
value of li^ and will provide little damping. 

Other important factors include the multiple coupled 
degrees of freedom of the AMD and the directional nature 
of the piezo damping material, both of which are covered 
in the follow section. 


B. Finite Element Model (FEM) 

The simple 1-D model introduced in the previous sec¬ 
tion is useful to provide an intuitive understanding of 


the AMD damping mechanism. However, it is inade¬ 
quate to represent the details of the interaction between 
the AMD and the test mass acoustic modes (TMAMs). 
Each AMD has at least six resonant modes and hundreds 
of TMAMs are present in the frequency band of interest; 
a Finite Element Model is required in order to properly 
reproduce these modes and calculate the mode overlap 
between them. 

A FEM of the Advanced LIGO test mass with two at¬ 
tached AMDs was constructed with the ANSYS program 
m- The AMD model corresponds to the parameters 
of our prototype AMD (see Table |^, and the test mass 
model parameters are reported in Table |VTj All dissipa¬ 
tion mechanisms in the test mass substrate, coating and 
bonds were included in the FEM, along with all the losses 
related to the acoustic mode damper structure (see table 
|V|for a full list). 

A piezoelectric material (PZT), for which energy dis¬ 
sipation can be easily controlled, is ideal for AMD con¬ 
struction. The AMD design modeled here has 2 PZTs 
sandwiched between a reaction mass and the interferom¬ 
eter test mass (see figure]^. The PZTs operate in shear. 
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TABLE L List of the components for acoustic mode damper prototype. The dimensions and loss angle value are used in finite 
element modeling described in Sec. |mB] 


Component 

Material 

Dimenssions 

Loss angle 

Resistor 

Ceramic, surface mount 

0.7x1.2x0.5 mm^ 

10 kQ 

Reaction mass 

Tungsten 

0.7x1.2x0.5 mm^ 

4e-C 

1 

Epoxy 

conductive TruDuct 2902 

25 fim 

O.IC 

1 

Piezo 

shear plate TRS 200HD 

0.7x1.2x0.5 mm^ 

0.014 

Base 

fused silica + gold coating 

0.7x1.2x0.5 mm^ 


£ 0.7' 

Epoxy 

non-conductive (EP30, MasterBond) 

25 /xm 

O.ljJ 

] 


^ Assumed loss value. 

^ Loss obtained from m- 



FIG. 3. A viscous-like behavior of the loss angle for shunted 
piezoelectric material. The total loss factor is a combination 
of the shunting loss and structural loss of the piezoelectric 
material. The maximum peak position is located at 40 kHz 
which corresponds to the shunt with 10 kQ resistor. The 
larger the resistor the lower is frequency of the peak loss. 
The peak height is proportional to electromechanical coupling 
coefficient ki^. Note, the loss angle at low frequency section is 
mainly dominated by the structural loss angle of piezoelectric 
material. 


and are poled orthogonally to ensure that all but 1 of the 
6 lowest frequency AMD modes are damped. See table 
for a list of modes; the compression mode is not damped 
by the shear plates. 

For the shear resonant damper, the spring constant can 
be associated with the shear deformation of the piezoma¬ 
terial 


kil + irj)=Re{cr,){l + iv^^^)^, (5) 

where C 55 is the shear stiffness matrix element, S is the 
area, and h is the height of the shear plate, respectively. 
The loss factor is induced by shunting the shear 
plate with a resistor. The active stiffness component C55 
becomes a complex quantity with a nonzero imaginary 
stiffness due to the shunt. The magnitude of 


is strictly related to the impedance of the shear plate- 
resistor circuit. The loss factor of the shunted shear plate 
can be defined as 


pzt ^ 

Re{cll) ■ 


(6) 


Because the impedance of any PZT is capacitive, the 
loss factor is frequency dependent. As such, careful 
selection of the shunting resistor and PZT dimensions 
can be used to maximize loss in the band of interest, 
as shown in figure A more detailed discussion of the 
piezoelectric loss angle can be found in Appendix [A| 

We validated the FEM of the AMD by computing its 
principle resonances, and comparing them with direct 
measurements performed on a prototype AMD. A total 
of five principle resonances were identified, with three 
types of modes: two flag, two anti-flag, and one rotation 
mode (see Table 0 . All these modes are characterized 
by large shear deformation for the double piezo configu¬ 
ration in the AMD. The sixth compression mode was not 
measured as it does not involve shear of the piezo plate. 
Table [11] shows the good agreement between the output 
of the model and measurements on the AMD prototype. 

Harmonic analysis (finite element analysis with an ex¬ 
citation at a fixed frequency) was conducted to estimate 
the Q-factor of each of the TMAMs. This approach al¬ 
lows us to include frequency dependent variables such as 
the shear plate stiffness 5 R(c 55 ), and its loss angle (see 
Appendix 0 . 

The modal Q-factor of each resonant mode of the sys¬ 
tem was calculated as 


Q(/n) 


Y.m^rn{fn)tCLn{(t)my 


(7) 


where fn is the frequency of the acoustic mode (see 
figure [^. The loss associated with each structural com¬ 
ponent is treated separately; is the modulus of the 
strain energy in the c omp onent and tan( 0 ^) is its 
loss factor (see tables]^ and VI). 

The modeled resonant frequencies of 12 TMAMs were 
then compared to the measured resonance frequencies 
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TABLE IL Verification of the finite element model for AMD. 
Five principle resonances (two flag, two anti-flag and one ro¬ 
tational mode) obtained with the model are compared to the 
measured values. The x, y suffix corresponds to the shear 
plate polarization direction. The principle resonances of the 
AMD were measured with a capacitive bridge circuit where 
one of the matching capacitors was the AMD prototype. The 
difference between flag and anti-flag modes corresponds to the 
location of the rotation axis about which reaction mass rocks. 
For the flag pole the rotation axis is at the bottom surface 
of the shear plate whereas for the anti-flag mode the rotation 
axis is at the height of mass center of the reaction mass. 

Mode type FEM [kHz] Measured [kHz] 


Fy-mode 

9.96 

9.77 

Fx-mode 

12.87 

12.61 

R-mode 

23.28 

24.13 

aFy-mode 

38.36 

37.39 

aFx-mode 

50.51 

48.86 


III), obtained according to methods described 
IVl The front face displacement of the modes 


(see table 

in section _ 

are shown in Fig. 

Note that the agreement between calculated and mea¬ 
sured Q-values increases with frequency. This may indi¬ 
cate an additional dissipation process which is missing in 
the model. Nevertheless, a small relative frequency offset 
A/ below 1% in Table III indicates a good agreement of 
the FEM with measured values. 


TABLE HI. The Q-factor of the test mass modes without 
AMDs installed. The 4th and 7th column correspond to the 
relative frequency noise Af% and Q-factor ratio between cal¬ 
culated and measured value, respectively. 


Mode # Freq. [Hz] 

FEM 

Freq. [Hz] 

Measured 

Af % 

Q-factor 

FEM Measured Ratio 

1 

8128.3 

8150.9 

0.3 

37M 

1.9M 

19.5 

2 

10391.1 

10418.1 

0.3 

63M 

14M 

4.5 

3 

12999.1 

12984.7 

0.1 

29M 

15M 

1.9 

4 

15101.4 

15047.2 

0.4 

56M 

17M 

3.3 

5 

15151.0 

15539.1 

2.5 

55M 

16M 

3.4 

6 

19487.0 

19544.9 

0.3 

30M 

7.0M 

4.3 

7 

20113.6 

20185.5 

0.3 

27M 

13M 

2.1 

8 

24824.0 

24901.8 

0.3 

32M 

16M 

2.0 

9 

26504.3 

26681.2 

0.7 

48M 

18M 

2.7 

10 

29767.4 

29699.5 

1.0 

18M 

15M 

1.1 

11 

30912.1 

31003.3 

0.3 

18M 

12M 

1.5 

12 

32664.1 

32743.2 

0.2 

14M 

13M 

1.1 



26681.(2) Hz 





32743.(2) Hz 


FIG. 4. Test mass drumhead modes for which Q-factor was 
measured. The color-code corresponds to the test mass front- 
face displacement amplitude. This figure was obtained from 
FEM analysis of the test mass model without AMDs. 


IV. EXPERIMENTAL RESULTS 


Several AMD prototypes were constructed, each con¬ 
sisting of six components (see Table |^: 

Reaction Mass: A 10 g tungsten cylinder which tunes 
the AMD principle resonances to frequencies above 10 
kHz is located on top of the shear plates. 

Shear Plate: Two piezoelectric shear plates, oriented 
with perpendicular polarizations are glued to the reaction 
mass and base with conductive epoxy Tra-Duct2902. The 
epoxy serves to electrically connect the PZT electrodes 
to the reaction mass and base. 

Base: The interface between flat shear plates and 
curved barrel of the test mass. The top flat surface is 
gold coated with separate sections to which shear plates 
are glued. The bottom surface is curved and match to 
the test mass barrel curvature. The base is made from 
fused silica and glued to the test mass with nonconduc- 
tive epoxy. 

Shunting circuit: One 10 kQ resistor for each shear 
plate is glued to the reaction mass with conductive 
epoxy. The circuit is closed with 100 /im diameter silver 
coated copper wire, which is soldered to the resistor and 
gold coated surface of the base. 

To measure Q-factors of the test mass modes, both 
with and without AMDs attached, a 16 m optical cavity 
was used, with the test mass forming the end mirror of 
the cavity. The optical cavity supports a 2 mm diameter 
Gaussian mode; the resonant beam probes the motion in 
the center of the test mass face, so that modes with an 
anti-node in the center of the mass are easily measured. 
A total of 12 modes were identified (Fig.[^ and measured, 
with the results shown in Table mi 

In order to measure the impact of the AMD on the 
Q of the TMAMs, several modes of the test mass were 
excited using electrostatic actuators m and their ring- 
down times observed. This measurement was repeated in 
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three configurations: with no AMD, with 2 AMDs, with 2 
AMDs which lacked resistive damping (shunt wires cut). 

Each acoustic mode was detected in the cavity locking 
error signal as a peak in the Fourier domain. After being 
excited, each mode amplitude was recorded as a function 
of time to estimate the decay time r. The Q-factor was 
determined from the decay time r, according to 


resonance is close to the TMAM frequency, equation 
indicates that the mode Q can be decreased by lowering 
the AMD loss. 

The large Q reduction for modes 7^3 and 7^4 is due 
the on-resonance interaction between TMAM and AMD, 
whereas for modes 7^5 and 7^12 the large Q reduction is 
due to the anti-node AMD location on the TMAM. 


Q = tt/ot. (8) 

where /o is the resonant frequency in Hz and r is the 
exponential decay time constant. 



FIG. 5. Two acoustic mode dampers (AMDs) attached to 
the barrel of a suspended aLIGO test mass. The AMDs are on 
opposite sides of the test mass, ~ 45° down from the midline. 


Measurements of the test mass mode Q-factor per¬ 
formed after installing two AMDs clearly show the sub¬ 
stantial damping capability of the AMD prototype, as 
reported in Table [Hi] (see also Table VII in Appendix for 
additional details). 

The results indicate that out of 12 modes, 11 are sup¬ 
pressed by at least factor of 10 and in some cases by more 
than two orders of magnitude. The relatively large dis¬ 
crepancy between model and measurement for mode 7^2 
can be explained by the AMD and TM interaction con¬ 
dition for this particular mode. The FEM predicts that 
mode 7^2 will be within 500Hz of the AMD Fy-mode at 
10 kHz, while the measured values give a separation of 
800Hz (mostly due to the AMD resonance being off), re¬ 
ducing the interaction between the AMD and TM mode. 

In the off-resonance interactions, which are more com¬ 
mon and set the lower limit to AMD performance, the 
discrepancy is generally less than a factor of a few. 

Surprisingly, modes 7^2 and 7^6 show a counterintu¬ 
itive behavior; opening the resistive circuit of the AMDs 
to decreases the TMAM Q-factor. However, since elec¬ 
trical circuit of AMD affects mechanical stiffness of the 
shear plates it is expected that principle resonances also 
changed when the circuit is opened. Also, if the AMD 


TABLE IV. Test mass mode suppression obtained with two 
AMDs. 


Mode 7 ^ Freq. [Hz] 

Damping 

factor 

Resistive 

contribution [%] 

1 

8150.9 

32.2 

12 

2 

10418.1 

31.8 

-57 

3 

12984.7 

441.2 

21 

4 

15047.2 

81.0 

9 

5 

15539.1 

>320 

>71 

6 

19544.9 

15.6 

-2 

7 

20185.5 

13.8 

41 

8 

24901.8 

5.2 

4 

9 

26681.2 

>360 

>0 

10 

29699.5 

23.8 

43 

11 

31003.3 

307.7 

68 

12 

32743.2 

>260 

>4 


V. ANALYSIS OF THE AMD THERMAL NOISE 

The AMD is designed to increase the mechanical 
damping of the test mass acoustic modes above 10 kHz. 
At the same time, the AMD must introduce minimal ad¬ 
ditional mechanical loss in the 10-1000 Hz band, where 
low mechanical loss is required to keep test mass ther¬ 
mal noise small [18]. Thus it is critical to calculate the 
thermal noise resulting from the AMDs in our overall 
evaluation of this PI mitigation technique. 

We used the FEM described above to calculate the 
thermal noise resulting from our experimental test of two 
AMDs attached to a test mass. The AMD thermal noise 
was calculated numerically at 100 Hz, the most sensi¬ 
tive part of the detection band, using Levin’s approach 
m- The energy dissipation per cycle was computed us¬ 
ing Eqn. |A15| for a pressure profile corresponding to the 
Advanced LIGO geometry (laser beam radius ujq = 5.5 
cm, incident on the front face of the TM). Results are 
shown in Table jVj 

The test mass thermal noise level of 5.2 x 10 ^^m/V^Hz 
corresponds to the design level for Advanced LIGO, and 
is dominated by optical coating loss [201121]. As Table iq 
shows, while our prototype AMD would contribute sig¬ 
nificantly to thermal noise it is not orders of magnitude 
above the more fundamental sources of thermal noise. 

Relative to the prototype device, our model points to 
several design and material improvements that can be 
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TABLE V. Thermal noise budget for test mass and acoustic 
mode damper prototype when attached to the test mass. The 
thermal noise was calculated for the laser beam spot size of 
(jj = 5.5 cm. 


Gomponent 

Thermal noise @ 100 Hz 

Test Mass: 


Substrate 

0.8 

Optical coating 

5.1 

Suspension ears 

0.0 

Ears bond 

0.6 

Total per TM 

5.2 

AMD: 


Reaction mass (RM) 

0.3 

Epoxy (RM-PZT) 

4.8 

Shear plate (PZT-X]^ 

^ 5.5 

Shear plate (PZT-Y^ 

' 5.5 

Epoxy Base-PZT 

7.2 

Base 

0.0 

Epoxy Base-TM 

6.0 

Total per AMD 

13 


^ A product of a structural+resistive loss angle. Note that a 
major contribution to that value comes from the structural loss, 
see Fig. 

made to significantly reduce the thermal noise impact. 
The major AMD thermal noise contributors are the epox¬ 
ies used to bond the AMD elements, and structural loss 
in the piezo shear plates. The former can be improved 
with lower loss epoxy and thinner bond layers. The lat¬ 
ter can be improved with a more suitable choice of piezo 
material. 

Other design elements can also be altered. The mass 
of the reaction mass can be reduced to lower the thermal 


noise without significantly affecting the acoustic mode 
damping performance. Another modification would be 
to avoid alignment of the piezo shear plate polarization 
with the laser beam axis, to minimize the contribution 
of resistive loss to the thermal noise. These and other 
design optimizations will be explored in a future paper. 


VI. CONCLUSION 

Acoustic mode dampers represent a simple yet effec¬ 
tive approach to damping parametric instabilities. The 
great advantage of this approach over active damping [9] 
is that many test mass acoustic modes are effected simul¬ 
taneously, and no further intervention is necessary. This 
is likely to be a critical feature in instruments that would 
otherwise suffer from multiple acoustic modes simultane¬ 
ously excited by parametric instabilities. 

The investigation presented here involved modeling 
and construction of a prototype AMD, which was shown 
to effectively damp test mass acoustic modes. The ther¬ 
mal noise associated with this prototype AMD was also 
computed. Though the prototype AMD does not meet 
the stringent thermal noise requirements of gravitational 
wave detectors, several design elements were identified 
for improvement, making this approach a viable solution 
to parametric instabilities. 
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Appendix A: Piezoelectric material 


Any piezoelectric material is strictly characterized by 
electromechanical properties. The fact that both electri¬ 
cal and mechanical properties are inseparable allows us 
to represent the dissipation process in a shunted piezo¬ 
electric material either as a Jonson heat or mechanical 
loss. For the purpose of this work we focus on disspation 
process in terms of the mechanical loss. 

From the stress-charge of the piezoelectric constitutive 
equation it is straightforward to derive the total induced 
stress in the shunted piezoelectric material [22] . Assum¬ 
ing no external current plugged to the piezoelement elec¬ 
trode we get 


ZTOT ^ ^. (A2) 

The admittance of piezelement is assumed to be ex¬ 
clusively capacitive. Hence, 

= iujAe^L~^ = iujC^ = iujC^ (e^) ^ (A3) 

where C, e is the capacitance and the dielectric constant 
matrices under constant strain S and constant stress T, 
respectively. It is convenient to operate with since 
this quantity can be easily measured at stress free condi¬ 
tions and no shortened piezelement electrodes. 
According to Eqn jAH if the piezoelement is integrated 
with a nonzero impedance electric circuit, the imaginary 
part stiffness tensor arises. The imaginary part can be 
interpreted as a dissipative component of the stiffness 
tensor. The total shunted stiffness matrix is 


= c^ -iuje^Z^^^AeL-\ 


(A4) 


Because matrices A, and L a re d iagonal and 

knowing that AL~^ = Eqn. A4 can be writ¬ 

ten as 


r,SU _ 


— e^eZ (e*^) , 


(A5) 


where Z is the nondemensional impedance 

Z = Z^^^Y^. (A6) 


Note, for zero impedance shunted circuit sets Z to be a 
unity matrix Z = I whereas for the nonzero impedance 
shunting Z is a complex quantity. Using indexing 
Eqn|A4| can be written as 


^SU — ^E _ ^ki^kj^kk 
^kk 


(A7) 


zTOT 


^(6x1) 


= c 


E 

(6x6) 


— lUJ- 


■^(6x3)^(3x3)^(3x3)e(3x6) 

L(3x3) 


^(6x1): 


(Al) 

where is the mechanical stiffness matrix under con¬ 
stant electric field, and e are the piezoelectric stress 
constant. The upper script t refers to transpose operator 
[23] . Electrode area is represented by matrix A, whereas 
thickness of the piezomaterial between electrodes by ma¬ 
trix L. The strain .S' is a product of the acting stress a 
on piezo-element and charge accumulation in the piezo¬ 
element. Each bracket corresponds to the matrix dimens- 
sion. It is assumed that piezoelectric element has a brick 
shape thus matrices Z^^^,A, and L are diagonal. 

The total impedance in Eqn |Al| is inverse sum of the 
piezoelement admittance under constant electric field Y^ 
and the admittance of the external circuit con- 


where k = 1,2,3 and corresponds to the electrode posi¬ 
tion, is the dielectric permittivity of the piezoelectric 
material. 

From Eqn |A3| and |A7| becomes clear that loss can be 
control with the shunting circuit. Moreover, such loss is 
in fact a frequency dependent quantity with arbitrary 
shape of the loss curve dependent on shunting circuit. 

For the shear plate with the single pair of electrodes 
only a stiffness matrix element C55 is affected by shunting 
resistor, thus 


„su _ „E 615615X11 

655 ~ 655 g 

^11 


(A8) 


and the nondemensional complex impedance Zkk has a 
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single nonunity element 


iujef^RC"^ 

^ iujef-^RC^ + efi ^ 

[u^ef.RCn" , 

[u^ef.RCTf + [el,f 


(A9) 


Since the stiffness matrix is a complex quantity we 
can define the loss factor r] as 


Re{c^^) ’ 


(AlO) 


and thus 


Vkij — 


(-^/c/c) X/c i jf 

Re{Zkk)Xkij + 1 


(All) 


where Xkij = ^ki^kj{cfje^^)~^. Index k is associated with 
the electrodes orientation, whereas indices i and j corre¬ 
spond to stress-strain directions in the stiffness matrix 
^su rj^Yie shunting loss factor of the shear plate becomes 


^(-^ll)Xl55 
5R(^1 i)X155 + 1’ 


(A12) 


This is the main loss mechanism based on which shear 
AMD operates. Note that according to Eqn |A8[ the 
material stiffness (the real part of is also reduced 
and should be included in analysis. 


1. Total mechanical loss angle 

An assumption of the stiffness matrix being real in the 
absence of shunt is not sufficient for accurate estimation 
of energy dissipation. It becomes especially important in 
the thermal noise analyisis, see Sec. |V| 

For known material loss angle of the nonshunted piezoele¬ 
ment, the total loss factor of piezoelectric material can 
be defined as 


where tan((/)^"^) is the material loss factor matrix of 
piezo-element. The loss factor tan(0^"^) of can be 
easily computed using Eqn. 

All piezoelectric material has anisotropic structure 
what implies that strain energy dissipation in such mate¬ 
rial must depend on the piezoelement geometric shape. It 
is more convenient to use effective loss angle 0e// rather 
than the loss angle 0 for each stiffness matrix compo¬ 
nent. We define the energy dissipation per cycle in the 
piezo-element 

Wdia = 2TrWatVeff (A14) 

= J Re{SiS*)Re{c^^)7]^fdV (A15) 

where Wst is the stored strain energy such that 2'KWst = 
/ Re{SiTj)dV, where S, T are the complex strain and 
stress, respectively and V is the volume of the piezo¬ 
element. 

Since loss factor of the piezoelement is inverse of its Q- 
factor or equally a ratio of dissipated energy per cycle 
Wdis to energy stored Wst in the piezoelement, we can 
defined the effective noise factor rj^ff 


_ Wdis _ f Re(SiS*}Re(c^jf)7jfddV 
27rWst jRe(SiT*)dV 


(A16) 


where S', T are the complex strain and stress, respectively 
and V is the volume of the piezo-element. This is a key 
equation in the finite element analysis which properly 
estimate contribution of piezoelemnt in the strain energy 
dissipation in the test mass. 

In our analysis we assumed constant material intrinsic 
loss factor for all elements equal to tan(0^"^)=O.O14 
[24] . what leads to 


0.014+ /m(Zii)xi55 
1 + Re(Zii)xi55 


(A17) 


Note, the remaining elements of shunting induced loss 
are equal to the material structural loss factor. For 
the shear plate configuration X155 is simply a function of 
electromechanical coupling coefficient ki^ and is equal to 
k^ 

Xi55 = i_k 2 > The total loss angle for the cfg stiffness 
element is shown in Figj^ 


^tot _ + Im{Zii)xkij 

~ 1 + Re{Zu)xkij 


(A13) 


Appendix B: Additional Tables 
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TABLE VL Values for the aLIGO end test mass parameters used in this paper. 


Optical parameters: 

End test mass transmissivity 5.0 ppm 

1 layer each of Ta 2 05 /Si 02 with thickness of 2.32 /rm, and 3.49 respectively 

Mechanical properties of optical coating and 

ear bond: 




Ta2 05 

Si02 

Ear bond 

Young’s modulus 

140 GPa 

70 GPa 

7.2 GPa 

Poisson ratio 

0.23 

0.17 

0.17 

Density 

8300 kg/m^ 

2201 kg/m^ 

2202 kg/m^ 

Refractive index 

2.06539 

1.45 

- 

Loss angl^ 

2.4-10-4+fT.8-10-9 

0.4-10-4+f-l.1-10-9 

0.1 

Test mass dimensions: 




Radius 

0.17 m 



Thickness 

0.2 m 



Flats width 

0.095 m 



Wedge angle 

0.07 deg 



Mass 

40 kg 



Loss angl^ 

7.6-10-i2.f 0.77 



Materials 

fused silica 




^ For the purpose of numerical analysis, a multilayer coating was reduced to double layers with a total thickness which corresponds to 
the sum of all 18/19 layers of Ta 205 /Si 02 for the ETM. 

^ f-frequency. Note, the difference in loss angles for the substrate, and for the optical coating made from fused silica. Losses were 
obtained from m- Additionally, coating loss angles were revised to the current measured values. 

^ Both test mass and suspension ears. 


TABLE VII. Results of the Q-factor measurement with at¬ 
tached AMDs for shunted and non-shunted cases. 


Q-factor 

Mode # Resistive Open Circuit 
EEM Measured EEM Measured 


1 

52k 

59k 

70k 

67k 

2 

7.9k 

440k 

9.9k 

280k 

3 

23k 

34k 

23k 

43k 

4 

420k 

210 k 

510k 

230k 

5 

2.4M 

<50k 

2.2M 

170k 

6 

1.9M 

450k 

4.8M 

440k 

7 

I.IM 

940k 

1.6M 

1.6M 

8 

6.7M 

3.1M 

9.5M 

3.2M 

9 

49k 

<50k 

64k 

<50k 

10 

1.9M 

630k 

3.3M 

I.IM 

11 

61k 

39k 

116k 

120 k 

12 

25k 

<50k 

260k 

52k 











